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Abstract 

D ■ Splines can be constructed by convolving the indicator function 

I of the Voronoi cell of a lattice. This paper presents simple crite- 

ria that imply that only a small subset of such spline families can 
be refined: essentially the well-known box splines and tensor-product 
splines. Among the many non-refinable constructions are hex-splines 
and their generalization to non- Cartesian lattices. An example shows 
how non-refinable splines can exhibit increased approximation error 
upon refinement of the lattice. 

1 Introduction 

Univariate B-splines are defined by repeated convolution, starting vi^ith the 
indicator functions of a partition of the real line by knots (An indicator 
function takes on the value one on the interval but is zero otherwise). This 
construction implies local support and a number of desirable properties (see 
[dBTHlldBHTj l that have made B-spl ines the representation of choice in mod- 
eling and analysis. In particular, B-splines can be exactly represented as a 
linear combination of B-splines with a finer knot sequence. This refinability 
is a key ingredient of multi-resolution, adaptive and sparse representation 
of data. 

By tensoring univariate B-splines, we can obtain on Cartesian grids in 
any dimension. For uniform knots, box-splines |dHR93| generalize this con- 
struction by allowing convolution directions other than the orthogonal ones 
of tensoring. This is not to say that box-spline convolution directions are 
arbitrary; to be practically useful, the directions need to be compatible with 
the lattice on which the spline is shifted, so that only a small number but 
sufficiently many lattice-shifts of the box-spline overlap at every point. 

As a most prominent example in two variables, the three direction box- 
spline forms a function with hexagonal footprint. The function is called 'hat 
function' and consists of six linear function pieces over the constituent trian- 
gles. Shifts over an equilateral triangulation add up to 1. Convolution of this 
hat function with itself results in a twice continuously differentiable function 
of degree 4 and m-fold convolution yields a function of degree 3m — 2 with 



1 



smoothness C^™'. Since this progression skips odd orders of smoothness, 
van der Ville et al. vBU"'"04] proposed to directly convolve the indicator 



function of the hexagon and build splines on the corresponding tessellation 
of the plane. They went on to show that the resulting hex-splines share a 
number of desirable properties familiar from box-splines. But the authors 
did not settle whether the splines were refinable [vdVUlO] . i.e. whether hex- 
splines of the given hexagonal tessellation T can be represented as linear 
combinations of hex-splines based on a finer-scale hexagonal tessellation, 
say ^T. Refinability is important in practice since it guarantees monotoni- 
cally decreasing approximation error as the scale of the tessellation refined. 
Moreover refinability is needed to locally adapt the space to features of 
higher frequency, a pre-requisite for multi-resolution analysis. 

• This paper presents simple criteria on a lattice that need to hold in 
order for shift-invariant functions on that lattice to be represented as 
linear combinations of piecewise constant shift-invariant functions on 
the smaller-scaled copy of the lattice. 

For example, the lattice must contain, for every of cell facet, the whole plane 
containing it. Therefore, requiring refinability, even of just the constant 
spline, strongly restricts allowable lattices. 

• In contrast to tensor-product and box splines, we show that hex-splines 
and similar constructions can only be scaled, but not refined: scaled 
hex-spline spaces are not nested. 

• A concrete example illustrates that non-refinable spaces can exhibit 
increased approximation error on a finer-scaled copy of the underlying 
tessellation (see Example [3]). 

• The analysis is extends to overcomplete spaces. 

Overview. Section [2] reviews lattices, hex-splines and their generaliza- 
tions. Section [3] exhibits two simple criteria for testing whether a lattice can 
support a refinable space of splines via convolution of indicator functions. 
Section m extends this investigation to overcomplete coverings by combining 
several families of indicator functions shifted by less than the lattice spacing. 
Section [5] illustrates the importance of refinability. 



2 Splines from lattice Voronoi cells 

A fc-dimensional lattice is a discrete subgroup of full rank in a /c-dimensional 
Euclidean vector space. A lattice may alternatively be viewed as a tessel- 
lation of space by identical cells. The Euclidean plane admits three highly 
symmetric tessellations into equilateral triangles, squares, or hexagons re- 
spectively. Convolution starting with the indicator function on any of these 
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polygons yields a hierarchy of spline functions of local support. The reg- 
ular partition into squares gives rise to uniform tensor-product B-splines 
and the regular triangulation and its hexagonal dual to box splines. An 
interesting and natural complement, to convolve the indicator function H 
of the hexagon with itself, was developed and analyzed by van De Ville et 



al. jvBU"'"04 . This yields a family of splines, of smoothness n — 1 sup- 
ported on a local n -|- 1-neighborhood, that the authors named hex-splines, 
van De Ville et al. observed that hexagons are Voronoi cells of a lattice 
and that the cell can be split into three quadrilaterals using one of two 
choices of the central split. Thus H can be viewed as the union of three 
constant box splines |dHR93] . an approach that was worked out more gen- 
erally for the FCC lattice by [KimOSa] and in more generality in [MElOj . 
[vBU+04 compares tensor-product splines and hex-splines and presents a 
Fourier transform. The transform allows determining an approxima- 
tion order with emphasis on low frequencies, as a combination of projection 
into the hex-spline space and a quasi-interpolation error. [CvBOSj derived 
quasi-interpolation formulas and showed promising results when applying 
hex-splines to reconstruction in image processing. 



3 Refinability constraints for lattice Voronoi cells 

Given a tessellation T of M'^, we denote by x{T) the space of indicator 
functions over the cells of T and by x{T^) the space of indicator functions on 
the smaller-scale copy of T. The space xi^) is refinable if each indicator 
function in x(^) can be represented as linear combinations of functions in 

Testing whether a tessellation admits a refinable space of functions re- 
quires off hand solving for weights such that a linear combination of elements 
in x{T^) with these weights reproduces each element in x{T). Proposi- 
tion [1] below provides a much simpler necessary condition that avoids such 
algebraic analysis. While we are interested in shift-invariant tessellations, 
Proposition [1] applies more generally. 

Proposition 1 LetT be a tessellation ofM.^ andT^ its smaller-scale copy. 
Then x{T) is refinable only if every tessellation facet of T is representable 
as a union of tessellation facets of . 

Proof Since is a tessellation, its cells do not overlap. Therefore, if a 
facet of a cell c in T is not a union of tessellation facets of then a cell 
c' of must cross this facet. Let G xi^^) be the indicator function of 
c' and H the indicator function of c. Then, in order to reproduce H, 
would have to take on both values and 1. ||| 

Scaling of the tessellation transforms this criterion to a less local one. 
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Proposition 2 For x{T) to he refinable, a tessellation T o/M'^ must con- 
tain, for each facet, a whole plane of the same dimension parallel to it. 

Proof Considering ever coarser tessellations, Proposition [T] implies that ar- 
bitrarily large extensions of each type of facet must be a union of tessellation 
facets of r. Ill 
The lattice structure implies that each such plane is replicated at all 
lattice points. 

Corollary 1 For x{T) to he refinahle, if T is a lattice, T must contain, for 
every facet, the whole plane of the same dimension that contains it. 

By inspection of the three regular tessellations of the plane, only the 
Cartesian grid and the uniform triangulation satisfy Corollary [H but not 
the partition into hexagons. 

Corollary 2 Hex splines are not refinable. 

We can generalize this observation by simplifying the inspection criterion. 




Figure 1: A pair of abutting facets, whose normals (dashed) have a strictly 
positive inner product, does not allow for refinability. 



Proposition 3 Consider a lattice T such that the reflection c' of a cell c 
across one of c's facets is again a cell ofT. If abutting cell facets of c meet 
with an obtuse angle then x(^) is not refinahle. 

Proof Let c be a cell of T with facets /i and /2 that join with an obtuse 
angle (see Fig. [1]). Let c' be the cell sharing /2. By the reflection symmetry, 
d also has an obtuse angle between /2 and the mirror image of /i across the 
plane F2 through /2. But on the side of c', Fi forms an accute angle with 
/2 and therefore intersects the interior of cell c'. Since cells can not be split, 
Fl cannot be part of T and the claim then follows from Corollary [TJ 1 1 1 
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Proposition [3] allows us to quickly decide which of the (regular crystallo- 
graphic) root lattices An, A^, Bn, T^n, V^, £j, j = 6, 7, 8 |CS98j are suitable 
for building splines by convolution of their Voronoi cells. 

Corollary 3 Splines obtained by convolution of Voronoi cells of regular 
crystallographic root lattices are not refinable, except for the Cartesian grid 
and the bivariate lattice with triangular Voronoi cells. 

Proof We test whether the Voronoi cells of the root lattices contain a pair of 
faces that meet with an obtuse angle. We may assume that one cell center 
is at the origin and take the inner product of the position vectors of two 
adjacent nearest neighbors, as identified by their root system; if the product 
is strictly positive, the corresponding Voronoi faces meet with an obtuse 
angle. 

The An lattice is traditionally defined via an embedding in R"^^, n > 1. 
Alternatively, Theorem 1 of [KPlOj gives a convenient geometric construc- 
tion in R" via the n x n generator matrix A„ := I„ + ^Jn, where I„ is the 
identity matrix, J„ the n x n matrix of ones and Cn ■= Denoting 
the ith coordinate vector by ej , we choose ei and ei + e2 on the Cartesian 
grid, and map them via A„ to the nearest An neighbors of the origin. The 
inner product of the images of ei and ei + e2 is 



A„ei • A„(ei + es) = l^Jli^ = -(n + - 1) > 0. 



.2 

n n 



For the An lattice, the computation is identical except that c„ := n"^^ • 

The inner product is ^^J^^^^ (n^ — 2n — 2 + 2-v/r2~+T)) > 0. 

For the T>n lattice, defined in n > 3 dimensions, the generator matrix is 



-Jn-l -1 



(see e.g. Section 7 of |KPllj ) and 

D„ei •D„(ei +e2) = 3 > 0. 

Since D~* is the generator of , the inner product for P* is 2. 

For Bn, the Cartesian cube lattice has an inner product identifying its 
spline constructions as potentially refinable (which indeed they are). On the 
other hand, splitting each the cube by adding a diagonal direction [Kim08b] 
yields the inner product ei • j = 1. 

For Sq, we select the root vectors (1,1,0,0,0,0) and (1, 1, 1, 1, 1, -\/3)/2 
with inner product 1. For Sj, we select the root vectors (1, 1, 0, 0, 0, 0, 0) and 
(1,1,1,1,1,1,^2)72 with inner product 1. For £§, we select the root vectors 
(1, 1,0, 0, 0, 0, 0, 0) and jsf^ with inner product 1. 

Equilateral triangular Voronoi cells in yield the inner product —-j-- 



5 



4 Overcomplete hex-spline spaces 



Since the evaluation of hex-splines by convolving three families of box splines 
already leads to a large number of terms, it is reasonable to investigate 
whether redundant superposition make hex-splines refinable. To test whether 
we can build refinable frames, let {Tj}j=i j be a family of tessellations ob- 
tained by shifting T by less than the lattice spacing so that their union 
covers M.^ J-fold. The next example makes this concrete for J = 3. 
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Figure 2: A lozenge-shaped pair of triangles (with markers and 1) is 
in the common support of three half-scaled translated copies of the grey 
hexagon and not in the support of other half-scaled hexagons. Since the 
pair straddles the boundary of the hexagon, any linear combination of the 
three indicator functions needs to be both 1 and 0. 

Example 1 Denote by a tessellation of the plane into unit-sized hexagons 
and by Ti and T2 its shifts by ^ri and (see Fig. 0j. Let H{x) be the 
indicator function of the unit hexagon of Ts centered at the origin. Consider 
the three indicator functions of hexagons of the ^-scaled tessellations shown 
in Fig. \^ Since the three hexagons supporting the three functions intersect 
in a pair of triangles, any linear combination of the functions has the same 
value on both triangles. But since the pair straddles the boundary of H{x) 
the value on one must be and the value on the other 1 implying that the 
joint space is not refinable. 

The example points to a simple extension of Corollary [TJ 

Corollary 4 Consider a family {Tj}j=i..j of tessellations each covering M'^. 
Then the space Uj=i jXi^j) is refinable only if the tessellation obtained 
by intersecting {Tj} does not contain a cell straddling a cell facet of any 
coarser-scaled copy of any Tj . 

The main argument concerning straddling cells applies to more general 
tessellations than shifts of a single tessellation T. 
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Extending the train of thought, the following Example [2] shows that 
families without straddling triangle pairs need not yield a refinable space of 
indicator functions either. 




Figure 3: Propagation of values via neighboring triangles that share an edge. 
This yields the contradiction that ao+bo+co = 1 (center) and 01+62+02 = 
(upper right) since the propagation implies Oj = oq, 6j = 60 and q = cq for 
all i. 



Example 2 Consider shifts 
H2{x) := H{x - 

of the indicator function H{x). The three corresponding tessellations now 
intersect only in single triangles so that the scenario of Corollary^ does not 
apply. However, an algebraic argument with a simple geometric interpreta- 
tion proves lack of refinability. 

We want to find scalar, real-valued coefficients Oi, bi and Ci, i = (^1,^2) S 
W? such that the following refinement equation holds: 

H{x)= aiH{2x-j,)+biH2{2x-ji) + CiH3{2x-ji), 




), H^ix) := H{x 



c := 



vr 

cos-, s :-- 
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sin- 
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We associate the coefficients with the center of its support hexagon. Observe 
then that, when two triangles share an edge and H{x) has the same value 
G 0, 1 on both triangles, then the coefficients at the two non-shared vertices 
must be equal. For example aQ+bQ+CQ = v = ai + b^+CQ implies oq = ai- As 
indicated by the arrows in Fig. 0, the coefficients are therefore propagated, 
separately inside and outside the support of H . This contradicts the refine- 
ment equation in that both ao + &o + co = 1 and 02 + 62 + ci = ao + bo + co = 
must hold. 

So even the two natural extensions to overcomplete spaces of shifted hex- 
sphnes do not afford refinabihty. 

The propagation argument generahzes to face-sharing simphces in any 
dimension. And it generahzes from binary to m-ary refinement. 



5 Importance of Refinability 

Why do we care about refinabihty and nestedness of spaces? Approximation 
order is weh-defined even for sequences of spaces that are not nested. For 
example, the elegant Fourier-based estimates of vBU"'"04] show that hex- 



splines resulting from m convolutions have, for low frequencies, an ap- 
proximation order of m. But approximation order is concerned with asymp- 
totic estimates. In practice one is more interested in predicting approxima- 
tion error. 

The following example shows why, for predicting the approximation er- 
ror, nested spaces are highly desirable. 

Example 3 Let T-L^ be the space of indicator functions over a regular tessel- 
lation by hexagons of diameter and such that, at each level of refinement, 
the origin is the center of one hexagon. Denote by H the indicator function 
in whose hexagon is centered at the origin. Let f be a function ob- 
tained by smoothing out H , say by a degree 3 Hermite interpolant, over a 
distance of at most 2e from the boundary of the hexagon. 

Then the L^ approximation error to f from T-L^ is approximately 2TTe, the 
integral over the smoothing region. However, since li} does not contain a 
linear combination that can replicate H (see Fig.^), the approximation 
error to f from is approximately 27r| >> 27re. That is the approximation 
error increases when refining the scale. 



6 Conclusion 

We identified simple necessary criteria for tessellations to admit a refin- 
able space of (convolutions of) indicator functions. Lattices, in particular, 
must contain, for every facet, the whole hyperplane containing it. With 
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(a) natural and man-made hex-tilings 



(b) half-scaled hex-tiling 



Figure 4: Examples for increased approximation error on a finer-scale 
tessellation, (a) Image credit to Simon Fraser and Fastfloors.com. (b) Non- 
nesting of the hex partition. 

Corollary [3] we observed that the increased isotropy of the Voronoi cells of 
non-Cartesian root lattices prevents refinability. Increased isotropy of the 
Voronoi cells is however the main reason for considering non-Cartesian lat- 
tices in the first place: they have higher packing densities leading to more 
efficient sampling [ PM62j . We observed that even overcomplete spaces ob- 
tained by natural superposition of shifted hexagonal tessellations fail to pro- 
vide refinable spaces from convolution of indicator functions. Finally, and 
here we omit the details, of the semi-regular tessellations of the plane, only 
3.6.3.6, the hex-tri-tessellation, satisfies the criteria of Proposition (TJ and 
while its indicator functions are refinable, generalizing the construction by 
convolution fails to yield a family of higher-order splines sharing all good 
properties of box-splines. 

In conclusion, if we want refinable classes of splines, remarkably few op- 
tions exist apart of box splines and tensor-product B-splines. Conversely, 
it should be noted, that adjusting and combining the families of symmetric 
box-splines on crystallographic root lattices, exhibited and analyzed for ex- 
ample jKPllj . does yield splines for any level of smoothness that obey the 
underlying symmetries. 

Acknowledgement Zhangjin Huang kindly worked out the first proof 
of non-refinability for Example [1] and Example [21 two scenarios I posed. The 
current versions shorten his arguments. 
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